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Abstract 



The Kelvin problem of an isotropic elastic space subject to a concentrated 
load is solved in a manner that exploits the problem's built-in symmetries so 
as to determine in the first place the unique balanced and compatible stress 
field. 
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1 Introduction 

This paper deals with a classical problem in linear elasticity, whose solution 
in terms of displacements was given by Lord Kelvin in a short paper dated 
1848 [l2j. The problem (see Figured]) consists in finding the equilibrium 



state of a linearly elastic, isotropic material body occupying the whole space 
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Figure 1: The Kelvin problem. 



and being subject to a point load. Kelvin's is the fundamental solution (in 
modern terminology, the Green function) later used by other XIX century 
elasticians to determine a number of strain nuclei (such as e.g. centers 
of dilatation and rotation), that is to say, of solutions to problems involv- 
ing unbounded elastic domains subject to force doublets (or dipoles) with or 
without moment, in various combinations, all inducing singular stress states. 
Singular stress states are also encountered when modeling the strongly lo- 
calized residual stress fields due to, say, gas trapped in small cavitites or the 
thermal mismatch associated with an inclusion. More recently, Kelvin solu- 
tion has been widely used as a basic ingredient in integral-equation methods 
for solving boundary- value problems in elastostatics (e.g., in boundary ele- 
ment methods [13]), 

Accounts of Kelvin problem, of different levels of clarity, completeness 
and purposes, are found in many expositions of the mathematical theory of 
classic elasticity: we recall the books by Love [6], Sokolnikoff |10| . Malvern 
[7], Benvenuto [2], Davis & Selvadurai [3], Barber p], and Sadd [9], but 
there are others. All these authors tackle the problem by solving the Navier 
equations for a displacement field having the representation constitutively 
implied by a tentative representation of the stress field in terms of a given 
scalar potential. We here propose a different approach, where the basic equa- 
tions of an elasticity problem — namely, the equilibrium, constitutive, and 
compatibility equations — are not combined preliminarily into one partial 
differential equation for displacement, as Navier did, but instead used se- 
quentially, in a semi-heuristic fashion, guided by the displacement and stress 
symmetries intrinsic to the problem at hand. 

In the two subsections to follow, we briefly sum up first the classical 
resolution technique, then ours; certain not completely obvious technical 
steps of the former, that are not found in any of the elasticity books quoted 
above, are contained in Appendix 1. Displacement and stress symmetries 
are discussed in Section 2; the balanced and compatible stress field that 
solves Kelvin problem is constructed in Section 3, the bulk of this paper; 
the easy task of deriving the associated strain and displacement fields is 
carried out in Section 4. 

1.1 Kelvin solution according to Love 

In his paper |12j . Kelvin does not disclose the details of the technique he 
adopted to solve the problem named after him: in little more than two pages, 
he just proposes the solution, stressing the similarity with the problem in 
the theory of thermal conduction, where a point heat source parallels the 
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role of a strain source under form of a concentrated load. We now give a 
quick account of Kelvin problem and its solution modeled after the account 
given by Love [6]. Love's procedure is essentially the same adopted by those 
authors who solved problems of the same type of Kelvin's after him; while 
the flow is clear, the mathematical developments are at times skipped. We 
warn the reader that the notation we use is rather different from the original 
- but coherent with the rest of the paper — and that some slight changes 
in presentation have been found either necessary or simply convenient. 
Navier equation for the displacement vector field u = u(x) is: 

(A + /i)V(Divit) + pAu + d = 0, (1) 

where A and [i are the Lame constants and d is the body force. With a view 
toward solving this equation, both the displacement and the force field are 
given a Helmhholtz representation in terms of potential pairs: 

u = V(p + curlw, d = V-0 + curl b, with Div w = Div b = oJl (2) 

The distance force is taken null in rl \ B p , where H is the whole space 
and Bp C rl denotes a sphere of radius p about the point o where the 
concentrated load / is applied. 

On recalling that AV(-) = VA(-) and Acurl(-) = curlA(-), equation 
(jl.ip can be written as: 

V((A + 2^)Av? + ^) + curl (/x Aw + VV> + b) = 0. (3) 

A solution of this equation can be obtained by solving the two equations: 

(\ + 2p)Atp + if; = 0, pAw + b = 0,E (4) 

2 We recall that, given any sufficiently smooth field u over a bounded regular region 1Z, 
its Helmhholtz representation consists in the pair of a scalar field ip and a divergenceless 
vector field w over 1Z such that 

u — V</3 + curl w with Div w — 0, 

(for u £ C(TZ) n C M (TZ), M > 1, both ip and w are of class C M (TZ)). A straightforward 
application of the identities curl X7tp = and divVy = Aip, yields: 

curl u — curl curl w and div u — Atp. 

3 Note that is recovered by taking the gradient of the first equation and the curl 

of the second, and then summing up the results. 
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both set over % \ B p ; the solutions ip p and w p of (jl.ip are: 
^ p (x) = ~ [ d(0-V € ( 7 - 1 (x,0)d«(0, 



where, for 7(x,£) := |x — £|, (47T7) -1 is the Green kernel of the laplacian. 
The potential pair yielding the representation ([Lip 1 of the solution to (jl.ip 
is found by taking the limits of (p p and w p for p — > 0: 

93 = — / -Vp, w = —f x Vp, p(x) := \x-o\. (5) 

Finally, the displacement vector is determined by inserting (jl.ip in (jl.iP j. 
For e2, e3 two unit vectors completing with the direction e\ of the applied 
load an orthonormal triplet, one finds: 

= (A+£)/ _/_ = / / 2(l-2i/) 1 x? 

Ul 8vr/z(A + 2/it) dx? A-npp 16vrG(l - 1/) \ p p p 3 

(X + p)f d 2 p fx x x 2 

U2 - 



U3 



S-Kp{\ + 2p)dxix 2 167rG(l-i/)' 

(X + p)f d 2 p fxixs 
8irp(X + 2p) dxix 3 16vrG(l - u) ' 

(6) 



where ?/,■:=«•&. 



1.2 Our method 

As shown in the previous subsection, it was Kelvin's concern to find a repre- 
sentation of the displacement field that would make Navier equations (jl.ip 
manageable; with his procedure, the elastic state {u,E,S} is determined 
sequentially, starting from the knowledge of u. We here propose a different, 
in fact inverse path. 

Firstly, we find the unique stress field S in % \ {o} that is: 

(i) balanced, in that it satisfies the balance equation 

DivS = 0, 

(ii) compatible, in that it satisfies the compatibility equation 

AS + — !— VV(trS) = 0, (7) 
1 + v 
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so that, in particular, 

A(trS) = 0; 



(8) 



(iii) balances the applied load, in the sense that 

/ Snda + f = 0, (9) 

for every full-dimensional neighborhood TV C H of the point o where 
the force / is applied (here n denotes the outer normal to the boundary 
dN of M). 

Once the stress field is found, the deformation field E can be computed 
by using the constitutive equation: 

E = h( s -ih {bS)I )' (10) 

where G and v are the shear and Poisson's moduli, and / is the identity 
tensor. Finally, the determination of the displacement field is accomplished 
by solving the differential equations: 

Vu + Vu T = 2E. 



2 Symmetries 

For a problem like Kelvin's, our intuition captures and restitutes a num- 
ber of symmetries that the displacement and contact-interaction fields must 
possess. As we shall see, recognizing such symmetries has three direct and 
useful consequences: 

(i) it reduces the representation of the stress tensor to a more tractable 
form; 

(ii) it makes it simpler to satisfy the compatibility equation; 

(iii) it furnishes a condition that is by all means equivalent to a bound- 
ary condition even in the present case of an infinite domain with no 
boundary. 

Different kinds of symmetries can manifest, depending on the body part we 
consider, either a right circular cylinder whose axis is directed as the load 
and passes through its application point (Figure [2]) or a ball centered at the 
load. 
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Figure 2: Cartesian axes and cylindrical coordinates. 



2.1 Symmetries of the displacement field 

We expect the displacement field that solves Kelvin problem to have cylin- 
drical symmetry, that is, when represented in the basis {e, h, h'}, to be such 
that both u z and u r are independent of <p and, in addition, is null: 

u(z,r,ip) = u z (z,r)e 1 + u r (z,r)h(ip). (11) 

Transforming this expectation into a parametric representation of the can- 
didate solution has two relevant consequences. 

1. Consequences on the strain and stress fields. In view of the definition 
of the strain measure E, we have that 

2 E^pz — i^tpix 'Uzhp ■ 

Accordingly, as it is not difficult to see, whenever (|2.ip prevail, both 
strain components E vr and E^ z turn out to be everywhere null. But 
then, given that "H is supposed to be comprised of an isotropic, linearly 
elastic material, both stress components S vr and Su, z must also be 
everywhere null: 

S vr = and S^ z = 0. 

Therefore, the class of stress fields of interest has the reduced repre- 
sentation: 

S = ei ei <g) ei + a 2 h ®h + a 3 h' <g) h' + a^{e\ (g> h + h <g) ei), (12) 
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parameterized by the four scalar- valued mappings (Xj, with 



0=1(2, 


r) 


— Szz — 


s ■ 


e x <B> ei , 




r) 




s ■ 


h <g> h, 




r) 


— S^ptp — 


■ s 


■ti®ti 


0=4(2:, 


r) 


= S zr = 


s ■ 


e\ ® /i. 



2. Consequences on the compatibility condition. It follows from (|2.ip that 
V« = u z , z e\®e\-\-u r , z h®e\-\-r~ u r h' '<g)/i' '+u z , r e\®h-\-u r , r h®h; 
consequently, 

E zz — U z , z , E rr — U r , r , E^tp — T U r , E zr — — ^U z , r ~\~U r ,z J ; 

(13) 

and hence, in particular, 

E rr = {rE w ), r . (14) 

In view of the inverse constitutive equation (|1.2j) , this consistency con- 
dition of representation (|2.1jl can be rewritten as: 

o"2 - (ro-3), r +vra, r = 0. (15) 

In addition to (|2.1|) . another symmetry condition prevails, namely, 

lim ti r (z, r) = 0; (16) 

as we shall see later on, this last condition works in all respects as a boundary 
conditions. 

Remark. Condition (j2J) is related to the compatibility equation 

curl curl E = 0. 

In fact, given the noted symmetries in the displacement and deformation 
fields, it is not difficult to see that the e\ <g> e\— component of curl curl E is 
given by: 

curl curl E ■ e\ ® e\ =: (curl curl E) zz = = r~ 1 E rr , r — 'E<pu>,rr ~ 2r _1 B™, in 
or rather, equivalently, 

a direct consequence of ([2]). 
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2.2 Symmetries on the traction 

Recall, to begin with, the following relationships between the cylindrical and 
spherical coordinates of a given point: 

z = pcost9, r = p | sin#|; p 2 = z 2 + r 2 , |tani?| = - 



(Figure[3]). 




Figure 3: Cartesian axes and spherical coordinates. 

When written as functions of the spherical coordinates, the stress com- 
ponents di will be denoted by <jj, with 

di(z,r) = d i (pcos'&,p\ sint?|) =: ai(p,$), 



a set of functions assumed to be such that 

*i(p,#) =3=i(p,-0) Vp>0. 



(17) 



Consider now a body part under form of a ball B p centered at the point of 
application of the load. On writing 

n(p, = costfei + | sini9| h(<p), (#,<p) € (0, +ir) x (0,27r) 

for the outer normal to the boundary dB p of B p , the traction on that surface 
has the expression: 

Sn = (cost?o"i + | sin$|a"4)ei + (| sin#|o"2 + cos & cr 4) h. 

According to (|mj) . the equilibrium condition is: 

„2 



fei 



Sn da 







2tt 







1 1 sin i? | d'd 



-2tt P 2 ( J (cos0cri(p,0) + |sini?|5 4 (/9,i?))|sini?|di?)ei, 



(18) 
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whatever p > 0. Therefore, for the right side to remain finite when p is 
chosen arbitrarily big, it is necessary that 

ai(p,0) = p- 2 Ti(0) + o(/r 2 ); 
a 4 (/>,0) = p- 2 r 4 W+o(p- 2 ). 

This result suggests the Ansatz: 

ai{p,-&) = p - 2 t0), n(#)=n(-#) (i = l,...,4), (19) 
which takes into account the parity conditions (|2.2p . 



3 The stress field 

Given that, when using cylindrical coordinates, 

Div S = (Sei), z +{Sh), r +r- 1 (Sh'), ip +r~ 1 Sh, 

a stress field in the class ([1]) is balanced for null distance forces if 

= r {aie 1 + a i h) u +{r(a 2 h + a 4 ei)), r +(a ?l h'), ip 
= (r<7i, z +(ro- 4 ), r )ei + (ra±, z +(rcr 2 ), r -a 3 )h, 

in T~L \ {0}. The choice of the parameter mappings is therefore restricted to 
those satisfying the following partial differential equations: 

rai, z +(r<7 4 ),r = 0, (20) 
rcn,z +(ra 2 ), r -cr 3 = 0. (21) 

Moreover, it is not difficult to see that, under the present circumstances, the 
vectorial compatibility condition (jn|) is equivalent to a system of four scalar 
equations, namely, 

Aax + a }ZZ = 0, 
Aa 2 - 2r~ 2 (o- 2 - cr 3 ) + a, rr = 0, . 
Aa 3 + 2r- 2 (a 2 -a 3 )+r- 1 a,r = 0, 1 ' 

A(T4 — r _2 cj4 + a, zr = 0, 

where 

a := (1 + u)~hT S = (1 + ^) _1 (cti + a 2 + ct 3 ) (23) 
must be harmonic: 

Aa = 0, (24) 
to satisfy (jn]). To solve ([3]), we propose the following sequential procedures: 
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1. to determine, a multiplicative constant apart, an appropriate solution 
a = a(z,r) of the Laplace equation ([3|); 

2. to integrate ([3]) 1 for 0*1, in the form: 

a 1 = -A~ 1 [a, zz ] + c 1 a, (25) 

where A -1 denotes the integral operator that formally inverts the 
laplacian, and c\ is a constant to be determined. 

3. to integrate ([3]) 4 for 04; 

4. to determine the fields ij 2 and (73, by solving the system of ([3]) 2 and 
©: 

cr 2 + cr 3 = A V,^]+c 2 a, 
0"2 - (rcr 3 ), r +vra, r = 0, 
where c 2 is a second constant, such that 

ci + c 2 = 1 + i/0 (27) 

Remark. For an alternative procedure, note the following consequence of ([3j) 
and ©: 

V(rcr 4 ) = -((rcr 2 ), r -03) ei - fr. (28) 

Were step 4. taken right after steps 1. and 2., the field 04 could be de- 
termined by integrating equation ([3]) along any regular curve C, arc- length 
parameterized, with tangent t, and going from a fixed point xq to the vari- 
able point x: 

(™-4)|* = - ((r<r 2 ),r -<T 3 )ei + mi )2 h) • i ds. (29) 

Remarkably, given the extreme points, the choice of the joining curve is 
irrelevant, due to a well known result in the theory of differential forms: for 
1Z open and star-shaped, and for u := oj\e\ + w 2 e 2 + 0^363 a vector field of 
class C 1 (7?.), the differential form uj = ui\dxi + uj2dx2 + co^dxs is exact if and 
only if curia; = in 1Z. In our case, the differential form and associated 
vector field under scrutiny are: 

w = ((m 2 ),r -0-3) dz + rai, z dr, u> = ((ra 2 ), r -03) ej + r0i, 2 h; 



'This condition is arrived at by adding (J3]) and ([3]^ and by taking into account of @. 
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consequently, 

curio; = (cri, zz -((rer 2 ), r -a 3 ), r )h'. 

Now, by differentiating Q with respect to z and ([3]) with respect to r, and 
by subtracting the resulting relations, we obtain: 

Cl,zz -((ro^jr -°"3),r = 0, 

which allows us to conclude that curio; = 0, and hence that cj is exact. 

3.1 Determination of a 

Relations ()2.2[) imply a preliminary representation for tr S : 

trS = p- 2 T($), t(0) = t(-0). 
For such a field to be harmonic: 

A(p- 2 r(#)) = 0, 

function r must satisfy the following ordinary differential equation: 
sin $ r" + cos $ r' + 2 sin $ r = 0, 

whose only regular solution is 

r($) = ro cos 

All in all, we are induced to choose 

tr S = To p~ 2 cos t?, (30) 

with ro a constant proportional to be applied load, to be determined later on 
(a method to arrive at (|3.ip without using (|2.2p is expounded in Appendix 
2). At this point, on recalling ([3]), we have: 

Q = Q0 p -2 COS With CKo = (1 + ^)" 1 T . 

3.2 Determination of o\ 

We start by showing how to solve the equation §3§i, that we here reproduce: 

Acri + a, zz = 0, 
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with the use of the information we gathered so far: that the unknown field 
has the preliminary representation ([272]) x : 

a 1 = a 1 (p,#)=p- 2 T 1 (#), n(tf) = ri(-tf), 

and that 

a = a(p, ■&) = aop~ 2 cos # . 

The equation we have to solve for f\ is: 

t" + cot § t[ + 2ti + 3a cos tf(2 cos 2 § - 3 sin 2 tf) = 0. (31) 



A particular solution of (|3.2p can be obtained by the method of parameter 
variation: 

t^VH^ocos 3 ??; 

moreover, the even solutions of the homogeneous equation associated with 
(1372]) are: 

r^(i?) = /3 cost?, /3 = a constant, 
a^tf) = p- 2 (^a cos 3 i? + /3 cos^). (32) 



in conclusion, 



3.3 Determination of a 4 

We have to solve ([3|) 4 for a field CT4 that as the representation (|2.2p i . namely, 

a 4 = Z 4 (p,#) = p- 2 T i ($), t 4 (0)=ti(-0), 
The function T4 is determined by the ordinary differential equation 

7-4 + cot #7-4 + (1 - cot 2 tf)r 4 + 3a |sintf|(4cos 2 t? - sin 2 tf) = 0, 
whose particular solutions are: 

Q 

rJV) = -a cos 2 tf|sini?|, (33) 
while the even solutions of the homogeneous associated equation are: 
74^(1?) = 7 I sini?|, 7 = a constant. 

We conclude that: 

3 

a±(p,<&) = p" 2 (-a cos 2 tf|sin#| + 7 |sintf|). (34) 
12 



As to the constants /5 and 7 , we note that the fields (|3.2p and (|3.3() 
satisfy the balance equation ([3]) only if 

Po = 7o5 

with this choice, the balance equation ([3]) is also satisfied. Furthermore, 
since the part-wise balance condition (|2.2p can now be written as 



f = -2vr / (costfri + | sin#|r 4 )| sinf?| dd, 
Jo 

we find that 

a + 2/3 = -^-. (35) 

Remark. Were we dealing with another classic problem in linear elasticity, 
Boussinesq's, the problem of a half space with a concentrated load perpen- 
dicular to the boundary, we could employ so far the same solving procedure. 
At this point, though, we would have had to satisfy the boundary condition, 
and this requires that /3 = 0; the reader is referred to [8] for details. 
Remark. On using cylindrical coordinates, we have that 

3 z 3 * 



2 ( z 2 + r 2)5/2 T ^0 ^ 2 + r 2)3/2 ' 

3 T T 

r ) = o a Q /..9 , 15^75 + Pott 



2 u ( z 2 +r 2)5/2 ^ u ( z 2 + r 2)3/2 ' 



a(z, r) = ao- 



'( z 2 +r 2)3/2' 

A consequence of the first and third of these relations is: 

3 z 3 (3 

= -«o , 2 , 2 , 5/2 + —a(z,r), 
2 (z z + r z ) / z ao 

Since 

3 z 3 

A" [a, zz ] = -^«o (z2 + r2)5/2 , (36) 
a comparison with ([3]) shows that 

ci = /3 /a . (37) 
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3.4 Determination of 02 and 03 

The functions 02 and u 3 we seek are to solve system ([3]), which, with the 
use of (|3|), (|3.3p . and (|3.3|) . can be rewritten as: 

3 z 3 z 
o"2 + 0-3 = — o«o— ^ + («o(l + f) - /3 )^, 

T ' Z 

o"2 - (ra 3 ), r -3a v-r = 0. 

P 5 

This system can be solved sequentially for g 3 and C2; it can be checked that 
solutions have the form: 

a 3 (z, r) = - (a (l - 2*,) - /? ) -J - (a (l - 2*) - 2/3 ) ^ + 

3 z 3 z 
£2(2, r) = --"o^r + («o(l + f ) - /3 )^r+ 

+ (a (l - 21/) -P )^+ (ao(l - 2*/) - 2/3 )^ " ^T' 
parameterized by an arbitrary function g(z)Jl 



3.5 Wrapping up 

At this point, we have satisfied the equilibrium and compatibility equations, 
and we have balanced the load. However, we have not yet completely found 
the solution, because the form of function g(z) and the value of the constant 
an are still wanted: we shall find this information by exploiting the symmetry 
property ([2.10 of the displacement field. 

By using the inverse constitutive law ([1.20 9. we find that 

*W = ((«o(l -»)-P )ji + Ml - 21/) - 2/?o)^3 " ^) , 

whence 

1 /* ZT Z q( z\ \ 

Ur = r£ w = -— (Ja (l - !/) - /3 )^3 + Ml " 2^) - 2/? ) ^3 " ■ 

(38) 

With problems formulated on a domain with no boundary, an incomplete, no 
matter if inessential, specification of the solution is to be expected, due to an 

5 Recall that the constant /3 can be determined by means of (|3.3[) in terms of ao and 
the magnitude of the applied load. 
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inevitable deficiency of conditions. This is not the case for the Kelvin prob- 
lem, where the deficit is covered by the symmetry condition (|2.ip . namely, 

lim u r (z, r) = 0, 

r->-0+ 

which plays the role of a Dirichlet boundary condition. Inserting (|3.5p in it, 
we deduce that: 

q (1 - 2v) - 2/3 = 0, g(z) = 0. 



The first of these relations, together with (|3.3p . yields the values of the 
constants «o and j3 : 

/ o /(i - M ^ 

ao = -Mi^7y ^ = -Mi^)' (39) 

We are now in a position to write the expression of the stress field: 

f ( z 3 z 
87r(l — v) \ p° p 6 

f ( ZT^ Z 

M*,r) = --—f r 3— + (1 - 2i/) 

87r(l — v) \ p° p 6 



/(I ~ 2i/) z 
8vr(l - v) p 3 ' 



(40) 



£3(2, r) 



Remark. In case the procedure delineated in Remark 2 is applied, the dis- 
covery of condition (|3.3I) is delayed. It is convenient to choose the curve C 
as the union of two different curves C = C\ U C2 ■ 

Ci := {x{z, r, ip) e H | z > z , r = r , <p = ip }, C 2 := {x(z, r,tp) € H\r > r ,ip = cp }, 

whose tangent vectors are, respectively, 

t x = ei, h = h(tp ). 
The integral in (|3|) can be determined as: 

u) ■ t = / {u} z e\ + U) r h) • t = / uj z (t,ro)dt + / 0J r (z,t)dt, 
JC Jc Jzo Jro 

with 

u z := (m 2 ),r -0-3 w r :=r(7 llZ . 
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It is not difficult to see that 

f f ( z 2 r 2 , s zr (1 - 2u)ri - 2(2 - v)zM 

u t = — -i 3^- - (1 - 2v)-r - ± '-2. — ^ ^-0_o 

Jc 8tt(1-i/) V P 5 V Pi 

where p : = p(zo,ro). On invoking equation ([3]), we have that 

^, . . / / z 2 r 2 , .zr {\-2v)ri- 2(2 -v)zlrl 
ra 4 (z,r)-r d 4 (z ,r ) = 3^- - (1 - 2v) 



87r(l — v) \ p 5 p 3 p\ 

whence, by choosing (zo,^o) = (0,0), we recover (|3.5p 1 . 



4 The strain and displacement fields 

To deduce the Kelvin strain field, we have recourse to the inverse constitutive 
equation (|1.2j) . We find: 



16vrG(l 


- v)p 5 ' 


/ 




16vrG(l - 




/ 


z 


16ttG(1 - 


v) p 3 ' 


/ 




16ttG(1 


- ^)p 5 ' 



(41) 



^- IMd-,)/ 2 ' 2 -*'^"-^''' 
As to the displacement field, given (13.5[) and (13. 5ft . we have that 

Wr = 167rG(l-z,)7 ; (42) 

moreover, on integrating the obvious combination of ([2]) 1 and (Jl|), we obtain 
the following preliminary representation of u z : 

f ( 2(\-2v) 1 \ 

U2 = 16ttG(1 - P 7 J + (r) ' 

With this, d2J) 4 , and (gj, we find that 

/i'(r) = h(r) = h ei, 

and we take ho = (cf. (|l.ip 1 ). 
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Remark. For ho ^ 0, the vector ho represents an arbitrary translation of 
the whole space in the vertical direction. Such an indeterminacy was to be 
expected, because the solution of an elasticity problem on a region without 
boundary is always determined to within a rigid displacement. In this spe- 
cific vertical translation of an arbitrary amount is the only type of 
rigid displacement compatible with the postulated symmetries. We regard 
disposing of such an arbitrariness by setting ho = as a covenient completion 
of the boundary conditions. 

Appendix 1. Developments complementing those 
in Subsection 11.11 

From the classic point of view, a concentrated force is the limit of a body 
force field having a shrinking support; according to the precise definition 
found in [11], a sequence {dn} of body force fields defined on an open neigh- 
borhood TZ of a point o tends to the load / concentrated at o if: 

(i) d, t e c 2 (ny, 

(ii) dn = su H\B rn (o), where {B Tn {o)} is a sequence of spheres of radius 
r n such that r n — > when n — > oo; 

(iii) lim / d n =/; 

(iv) the sequence { J~ \ dn\} is bounded. 

Given the vector d, the determination of the fields ip and b satisfying the 
Helhmoltz's decomposition (jl.lP a can be achieved in two steps. 

1. On applying the divergence operator to the equation (jl.ip . we obtain 
the Poisson equation 

AV> = Div d, 

whose solution has the well-known representation (see e.g. [4]): 
■4>{x) = - I G(x,f)Divd(Od«(f)» x e H, 

where 

G(x,0 = (4vrp(x,e))" 1 , P (x,0 := \x - £| . 
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On recalling the identity 

Div (ipv) = ipDiv v + v ■ Vip, 

for (p and v smooth scalar-valued and vector-valued fields, and on 
using the divergence theorem, we obtain: 

J j3 r 

2. On taking the curl of d, we obtain: 

curl d = curl curl b = V(Div 6) — A 6 = —Ab. 
Thus, we have to solve another Poisson equation: 

— A6 = curl d, 

whose solution is: 

b(x) = ±- / G(x,t)coxld(t)dv(0. 

An application of Stokes theorem combined with the identity 
curl {ipv) = (pcwxl v + V99 x v 

yields: 

br(x) = -^J B d (0 x V^p-\x,0)dv(C)- 

We now compute the limits of ip r and b r for r — > 0, under the assumption 
that 



lim / d(x)dv(x) 



We find: 

ip(x) = —^—f ■ V(r~ 1 ), b(x) = — — f x V(r _1 ), r(x) := \x — o\ , 

47T 47T 

or rather, since 2V(r~ 1 ) = AVr, 

ijj(x) = ~A(f ■ Vr), b(x) = ~^{f x Vr). 
We can now write system as follows: 

A {»-l f w ) =0 ' 

A [ w - —f x Vr ) = 0, 

V ^ J 

and read out the particular solution (jl.ip . 
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Appendix 2. An alternative way to obtain ( 13.11 ) 

A way to find solutions independent of p to the Laplace equation in cylin- 
drical coordinates: 

Aa = a, zz +a, rr +r~ 1 a, r = 
consists in looking for solutions, if any, having the form 

a(z, r) = p a z b r c . 

It is not dfficult to see that, for such an Ansatz to be successful, the expo- 
nents a, b, c must satisfy the following condition: 

a(a + 26 + 2c + 1) + b(b - l)p 2 z~ 2 + c 2 p 2 r~ 2 = Vz, r > 0, 

or rather, equivalently, they must be chosen so as to satisfy the following 
three algebraic conditions: 

c = 0, 6(6-l) = 0, a{a + 2b + 2c+ 1) = 0. 

There are four possibilities: (i) b = and a = 0; (ii) 6 = and a = 1; (iii) 
6 = 1 and a = 0; (iv) 6=1 and a = —3; accordingly, the desired field must 
have the following form: 

~ z 1 

a(z,r) = Qo^r + a.\- + a 2 z + a 3 . (43) 

This result can be applied when, as it happens when dealing with Kelvin 
problem, the trace of a stress field that solves Kelvin problem is found 
to be proportional to a harmonic field that, when expressed in cylindrical 
coordinates, does not depend on p. In such an application, one starts from 
the representation ^ and quickly sets both constants a 2 and to zero, so 
as to comply with the physical palusibility requirement that the stress field 
- and hence its trace - vanishes at infinity. Next, one does the same for a±, 
this time on the basis of an application of a result due to A. Signorini, that 
we now recall in a version appropriate to our present context (cf. e.g. Sect. 
18 of 0). 

Signorini's Lemma. Let S be a stress field that balances the distance and 
contact forces d and c acting on a domain 1Z with boundary dlZ: 

div S + d = in 1Z, Sn = c on dTZ. 
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Moreover, let w be a conveniently smooth vector field on 1Z U dlZ. Then, 

/ (Vw)Sdv = / w(g>ddv+ w®cda. (44) 
Jn Jn JdK 

We specialize (jH) for 1Z = B p , a ball of radius /) centered at the point of 
application of the load, d = 0, and w = x; by taking the trace of the 
resulting identity, we obtain: 

/ tr S dv = / pn ■ Sn da. (45) 

Now, 

/ pSn ■ nda = O(p), 
JdBp 

because we have from (|m|) that 



/ Snda = 0(1). 

JdBn 



'dB p 

On the other hand, as to the left side of (jU), we have that 
r rp 

/ tr S dv = 27r(l + u) / / (ao cos i9 + sai) d-d ds; 
JBp Jo Jo 

for it to be O(p) as well, a\ must be set equal to 0. 
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